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1 Introduction 

The purpose of this paper is to give a variational multisymplectic formu- 
lation of continuum mechanics from a point of view that will facilitate the 
development of a corresponding discrete theory, as in the PDE Veselov for- 
mulation due to Marsden, Patrick, and Shkoller [1998]. This discrete theory 
and its relation to finite element methods will be the subject of a companion 
paper, Marsden, Pekarsky, Shkoller, and West [2000]. 

In this paper, we restrict our attention to non-relativistic theories, but 
on general Riemannian manifolds. The relativistic case was considered in 
Kijowski and Tulczyjew [1979], where the authors take an alternative ap- 
proach of inverse fields, effectively exchanging the base and fiber spaces. See 
also Fernandez, Garcia and Rodrigo [1999] .Q 

1 There are a number of reasons, both functional analytic and geometric for motivating 
a formulation in terms of direct particle placement fields rather than on inverse fields. 
For example, in the infinite dimensional context, this is the setting in which one has 
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Two main applications of our theory are considered — fluid dynamics and 
elasticity — each specified by a particular choice of the Lagrangian density. 
The resulting Euler-Lagrange equations can be written in a well-known form 
by introducing the pressure function P and the Piola-Kirchhoff stress tensor 
V (equations ( 2.18Q and (|2.21j ) below, respectively). 



We only consider ideal, that is nonviscous, fluid dynamics in this paper, 
both compressible and incompressible cases. In the compressible case, we 
work out the details for barotropic fluids for which the stored energy is a 
function of the density. These results can be trivially extended to isentropic 
(compressible) fluids, when the stored energy depends also on the entropy. 
Both the density and the entropy are assumed to be some given functions in 
material representation, so that our formalism naturally includes inhomoge- 
neous ideal fluids. However, in our discussion of symmetries and correspond- 
ing conservation laws considered in §|| we restrict ourselves, for simplicity 
only, to fluids that are homogeneous in the reference configuration. We 
elaborate on this point below. 

For the theory of elasticity we restrict our attention to hyperelastic ma- 
terials, that is to materials whose constitutive law is derived from a stored 
energy function. Similarly, we assume that the material density is some given 
function which describes a possibly heterogeneous hyperelastic material. 

A general formalism for treating constrained multisymplectic theories is 
developed in §H[ Often, constraints that are treated in the multisymplectic 
context are dynamically invariant, as with the constraint divE = in elec- 
tromagnetism (see, for example, Gotay, Isenberg, and Marsden [1997]), or 
divE = p for electromagnetism interacting with charged matter. Our main 
example of a constraint in this paper is the incompressibility constraint in 
fluids, which, when viewed in the standard Eulerian, or spatial view of fluid 
mechanics is often considered to be a nonlocal constraint (because the pres- 
sure is determined by an elliptic equation and, correspondingly, the sound 
speed is infinite), so it is interesting how it is handled in the multisymplectic 
context, which is, by nature, a local formalism. 

In the current work we restrict our attention to first-order theories, in 
which both the Lagrangian and the constraints depend only on first deriva- 
tives of the fields. Moreover, we assume that time derivatives do not enter 
the constraints, which corresponds, using a chosen space-time splitting, to 

the deeper geometric and analytical properties of the Euler equations and related field 
theories, as in Arnold [1966], Ebin and Marsden [1970], Shkoller [1998], and Marsden, 
Ratiu and Shkoller [2000]. Moreover, the relativistic approach adopted in Kijowski and 
Tulczyjew [1979] cannot describe an incompressible fluid or elasticity because the notion 
of incompressibility is not defined in the relativistic context. 
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holonomic constraints on the corresponding infinite-dimensional configura- 
tion manifold in material representation. We briefly discuss the issues re- 
lated to extending this approach to non-holonomic constraints and to space- 
time covariant field theories in the last section. 

Symmetries and corresponding momentum maps and conservation the- 
orems are considered in a separate section (§||) since they are very different 
for different models of a continuous media, e.g. homogeneous fluid dynamics 
has a huge symmetry, namely the particle relabeling symmetry, while stan- 
dard elasticity (usually assumed to be inhomogeneous) has much smaller 
symmetry groups, such as rotations and translations in the Euclidean case. 
We emphasize that although the rest of the paper describes general hetero- 
geneous continuous media, the results of §54. only apply to fluid dynamics 
that is homogeneous in the reference configuration (e.g., the fluid starts out, 
but need not remain, homogeneous), where the symmetry group is the full 
group of volume-preserving diffeomorphisms T>^. However, these results can 
be generalized to inhomogeneous fluids, in which case the symmetry group 
is a subgroup T>^ C T>^ that preserves the level sets of the material density 
for barotropic fluids, or a subgroup D^' cn c T>^ that preserves the level sets 
of the material density and entropy for isentropic fluids. This puts us in the 
realm of a multisymplectic version of the Euler-Poincare theory - one needs 
to introduce additional advected quantities as basic fields to handle this sit- 
uation (see the discussion on symmetry and reduction in ^). We remark 
also that all continuum mechanics models should satisfy material frame in- 
difference principle, which, as is well known, can be readily accomplished by 
requiring the stored energy function to be a function of the Cauchy-Green 
tensor alone (see, e.g. Marsden and Hughes [1983], Lu and Papadopoulos 
[1999]). 

We finally remark on the notation. The reader is probably aware that 
typical fluids and elasticity literatures adhere to completely different sets of 
notations, which both differ substantially from those adopted in multisym- 
plectic theories. In our notations we follow Gotay, Isenberg, and Marsden 
[1997]. The companion paper Marsden, Pekarsky, Shkoller, and West [2000] 
uses primarily notation from Marsden and Hughes [1983] and concentrates 
on models of continuum mechanics in Euclidean spaces and their variational 
discretizations. 
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2 Compressible Continuum Mechanics 



To describe the multisymplectic framework of continuum mechanics, we 
must first specify the covariant configuration and phase spaces. Once we 
obtain a better understanding of the geometry of these manifolds we can 
consider the dynamics determined by a particular covariant Lagrangian. 

2.1 Configuration and Phase Spaces 

The Jet Bundle. We shall set up a formalism in which a continuous 
medium is described using sections of a fiber bundle Y over X; here X is 
the base manifold and Y consists of fibers Y x at each point x € X. Sections of 
the bundle itxy :Y—>X represent configurations, e.g. particle placement 
fields or deformations. 

Let (B,G) be a smooth n-dimensional compact oriented Riemannian 
manifold with a smooth boundary and let (M, g) be a smooth iV-dimensional 
compact oriented Riemannian manifold. For the non-relativistic case, the 
base manifold can be chosen to be a spacetime manifold represented by the 
product X = B x K of the manifold B together with time; (x, t) S X. Let 
us set x° = t, so that x M = (x l , x°) = (x l , t), with a = 0, . . . ,n, i = I, . . . ,n, 
denote coordinates on the (n + l)-dimensional manifold X. Construct a 
trivial bundle Y over X with M being a fiber at each point; that is, Y = 
X x M 3 (x, t, y) with y € M — the fiber coordinate. This bundle, 

ttxy ■ Y -» X; (x, t, y) i-» (x, t), 

with ttxy being the projection on the first factor, is the covariant configura- 
tion manifold for our theory. Let C = C°°(Y) be the set of smooth sections 
of Y. Then, a section <j> of C represents a time dependent configuration. 

Let y a , i = 1, . . . , N denote fiber coordinates so that a section eft has a 
coordinate representation (p(x) = (x^, <p a (x)) = y a ). The first jet bundle 
J l Y is the affine bundle over Y whose fiber above y € Y x consists of those 
linear maps 7 : T X X — > T y Y satisfying Tttxy 07 = Id,T x x- Coordinates 
on are denoted 7 = y a , v a ^). For a section cp, its tangent map at 
x £ X, denoted T x <f>, is an element of J Y^r x y Thus, the map x 1— > T x cj) is a 
local section of J^Y regarded as a bundle over X. This section is denoted 
j 1 ^ and is called the first jet extension of </>. In coordinates, is given by 
(x M , (j) a (x),d^<j) a ), where do4> a = d t <j) a and dk<fi a = d x k(j) a . 

Notice that we have introduced two different Riemannian structures on 
the configuration bundle. The internal metric on the spatial part B of the 
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base manifold X is denoted by G and the fiber, or field, metric on M is 
denoted by g. There are two main cases, which we consider in this paper: 

(i) fluid dynamics on a fixed background with fixed boundaries, when 
B and M are the same and the fiber metric g coincides with the 
base metric G; a special case of this is fluid dynamics on a region in 
Euclidean space; 

(ii) elasticity on a fixed background, when the metric spaces (B,G) and 
(M,g) are essentially different. 

Both approaches result in background theories. The case of relativistic fluid 
and elasticity was considered by Kijowski (see, e.g. Kijowski and Tulczyjew 



We shall see later that its pull-back by a section <j) has the interpretation of 
the Jacobian of the linear transformation D<j>t . 

A very important remark here is that even though in fluid dynamics 
metrics g and G coincide, i.e. on each fiber Y x , g is a copy of G, there is 
no cancellation because the metric tensors are evaluated at different points. 
For instance, in (|2.1[) g(y) does not coincide with G(x) unless y = x or both 
metrics are constant. Hence, only for fluid dynamics in Euclidean spaces, 
can one trivially raise and lower indices and drop all metric determinants 
and derivatives in the expressions in the next sections. 

The Dual Jet Bundle. Recall that the dual jet bundle J l Y* is an affine 
bundle over Y whose fiber at y £ Y x is the set of affine maps from J l Y to 
A n+ X x , where A n+l X denotes the bundle of (n+ l)-forms on X. A smooth 
section of J l Y* is an affine bundle map of J l Y to K n+1 X covering ~kxy- 
Fiber coordinates on J l Y* are (n,p a ^), which correspond to the affine map 
given in coordinates by v a ^ t— > (n + p a ^v a ^)d n+l x. 

To define canonical forms on J 1 !^*, another description of the dual bun- 
dle is convenient. Let A = A n+1 Y denote the bundle of (n + l)-forms on Y, 
with fiber over denoted by A y and with projection ttya '■ A — > Y. Let 

Z be its "vertically invariant" subbundle whose fiber is given by 





(2.1) 



Z y = {z £ Ay | v J w J z = for all v, w 6 V y Y} 



G 



where V y Y = {v £ T y Y \ Tttxy • v = 0} is a vertical subbundle. Elements 
of Z can be written uniquely as 

z = Ild n+1 x + p a ^dy a A d n Xfl 

where d n x fM = d^d n+1 x, so that (x fl ,y a ,Il,p a tl ) give coordinates on Z. 

Equating the coordinates (x^, y a , II, p/) of Z and of J X Y* defines a vec- 
tor bundle isomorphism Z <-> J l Y* . This isomorphism can also be defined 
intrinsically (see Gotay, Isenberg, and Marsden [1997]). 

Define the canonical (n + l)-form ©a on A by @\(z) = (ity\z), where 
z 6 A. The canonical (n + 2)-form is given by Oa = —dQ\. If zaz : Z — ► A 
denotes the inclusion, the corresponding canonical forms on Z are given by 
= «az® a an d ^ = — rf@ = ^az^a- ^ n coordinates they have the following 
representation 

9 = p a ^dy a A d n Xfl + Ud n+1 x, n = dy a A dp/ A (f% — (ill A d n+1 x. 
Ideal Fluid 

We now recall the classical material and spatial descriptions of ideal (i.e., 
nonviscous) fluids moving in a fixed region, i.e., with fixed boundary condi- 
tions. We set B = M and call it the reference fluid container. A fluid 
flow is given by a family of diffeomorphisms rjt : M — > M with t/q = Id, 
where ryt(M) is the fluid configuration at some later time t. Let x € M 
denote the original position of a fluid particle, then y = 7ft(x) E M is its 
position at time t; x and y are called material and spatial points, respec- 
tively. The material velocity is defined by V(x,t) = (d/dt)r/t (x). The 
same velocity viewed as a function of (y, t) is called the spatial velocity. 
It is denoted by u; that is, u(y,t) = V(x(y),t), where x = (y), so that 
u = V o^ 1 =l)o))- 1 , 

Thus, in the bundle picture above, the spatial part of the base manifold 
B C X has the interpretation of the reference configuration, while an extra 
dimension of X corresponds to the time evolution. All later configurations 
of the fluid are captured by a section <f> of the bundle Y, which gets the 
interpretation of a particle placement field. Pointwise this implies that x in 
the base point (x, t) represents the material point, while y 6 Y( x ,t) represents 
the spatial point and corresponds to a position y = <p(x,t) = r] t (x) of the 
fluid particle x at time t. 
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Elasticity 



For the theory of elasticity (as well as for fluids with a free boundary), the 
base and fiber spaces are generally different; (B,G) is traditionally called 
the reference configuration, while (M, g) denotes the ambient space. For 
classical 2 or 3 dimensional elasticity, M and B have the same dimension, 
while for rods and shells models the dimension of the reference configuration 
B is less than that of the ambient space. 

For a fixed time t, sections of the bundle Y, denoted by <j)t, play the role 
of deformations; they map reference configuration B into spatial configu- 
ration M. Upon restriction to the space of first jets, the fiber coordinates v 
of 7 = (x,y,v) £ J^Y become partial derivatives d(f) a /dx^; they consist of 
the time derivative of the deformation <j> a and the deformation gradient, 
F a i = d(p a jdx % . The first jet of a section <f> then has the following local 
representation = ((x, t), 4>(x, t),(f)(x, t), F(x, t)) : X — > J 1 !". 

Using the map 0, one can pull back and push forward metrics on the 
base and fiber manifolds. In particular, a pull-back of the field metric g on 
M to B C X defines the Green deformation tensor (also called the right 
Cauchy-Green tensor) C by C* = (f>t(g), while a push- forward of the base 
metric G on B C X to M defines the inverse of the Finger deformation 
tensor b (also called the left Cauchy-Green tensor): c = ft -1 = ((pt)*(G). In 
coordinates, 



Cij(x,t) = g ab F a i F i j (x ) t), c ab (y) = Gij(F~ 1 ) i a (F~ 1 y b (y), (2.2) 



where F is thought of as a function of y. We remark that C is defined 
whether or not the deformation is regular, while b and c rely on the regularity 
of (j>t- Another important remark is that operations flat b and sharp (J are 
taken with respect to the corresponding metrics on the space, so that, e.g 



Notice that J restricted to the first jets of sections is the Jacobian of 
D(f>t, that is, the determinant of the linear transformation Dcftt', it is given 
in coordinates by 



It is a scalar function of x and t, invariant under coordinate transforma- 
tions. Notice also that J(x, t) > for regular deformations with 4>(x, 0) = 
x, F(x, 0) = Id because J(x, 0) = 1. 
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2.2 Lagrangian Dynamics. 



To obtain the Euler-Lagrange equations for a particular model of a contin- 
uous medium, one needs to specify a Lagrangian density C. Naturally, it 
should contain terms corresponding to the kinetic energy and to the po- 
tential energy of the medium. Such terms depend on material properties 
such as mass density p as well as on the constitutive relation. The latter is 
determined by the form of the potential energy of the material. We remark 
that such an approach excludes from our consideration non-hyperelastic ma- 
terials whose constitutive laws cannot be obtained from a potential energy 
function. 



Lagrangian Density. Let the mass density p : B — > M be given for a par- 
ticular model of continuum mechanics. The Lagrangian density C : J l Y — > 
A n+1 X for a multisymplectic model of continuum mechanics is a smooth 
bundle map 



£(7) = L( 7 )d n+1 x = K — P = - ^/fet\G]p(x) 9ab v a v b d n+1 x 

- y/d^\G]p(x)W(x, G(x),g(y),v a j )d n+1 x, (2.3) 

where 7 6 and W is the stored energy function. The first term 

in fl2.3|) corresponds to the kinetic energy of the matter when restricted to 
first jet extensions, as v a o becomes the time derivative dt4> a of the section 
<j). The second term reflects the potential energy and depends on the spatial 
derivatives of the fields (upon restriction to first jet extensions), i.e. on the 
deformation gradient F. 

A choice of the stored energy energy function specifies a particular model 
of a continuous medium. While different general functional forms distinguish 
various broad classes of materials (elastic, fluid, etc.), the specific functional 
forms determine specific materials. Typically, for elasticity, W depends on 
the field's partial derivatives through the (Green) deformation tensor C, 
while for Newtonian fluid dynamics, W is only a function of the Jacobian J 

(B 



Legendre Transformations. The Lagrangian density (2J3) determines 
the Legendre transformation ¥C : J l Y — > J x y*. The conjugate momenta 
are given by the following expressions 

Pa = -^ = pg ab v b oVd^G], pJ = ^ = - P ^Vd^G], (2.4) 
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n = l 



dL a 

dv% v " 



1 a b Trr dW 

- 2 9a^ov\-W + — 



Define the energy density e by 



Vdet[G] 



e = p a °v a o — L or, equivalently ed n+L x = K + 



(2.5) 



then 



n 



y/det[G\t 



The Cartan Form. Using the Legendre transformation ( |2.4[ ), we can pull- 
back the canonical (n + l)-form from the dual bundle. The resulting form 
on J^Y is called the Cartan form and is given by 

k / dW , 

Qc = pg a bV b oVtet\G]dy a A d n x - p-— v ^t\G]dy a A d n Xj 



+ 

-dS c . 



dv a j 
1 , 



dv a i V 3 



Vdet[G](T +1 x. (2.6) 



We set f2,c 

Theorem 2.1 below provides a nicer method for obtaining the Cartan 
form via the Calculus of Variations and remains entirely on the Lagrangian 
bundle J l Y . Moreover, the variational approach is essential for the Veselov 
type discretization of our multisymplectic theory. We present it here for 
the benefit of the reader, but remark that it is not essential for our current 
exposition and can be omitted on a first reading (see Marsden, Patrick, and 
Shkoller [1998] for details). 



Variational Approach. To make the variational derivation of the equa- 
tions of motion rigorous as well as that of the geometric objects, such as the 
multisymplectic form and the Noether current, we need to introduce some 
new notations (see Marsden, Patrick, and Shkoller [1998]). These are gener- 
alizations of the notations used in the rest of the paper. They only apply to 
the variational derivation described here and later in Subsection [O] and do 
not influence the formalism and results in the rest of the paper. The reason 
for such generalizations is very important yet subtle: one should allow for 
arbitrary and not only vertical variations of the sections. 

Vertical variations are confined to the vertical subbundle VY C TY, 
V y Y = {V € TyY\TirxY • V = 0}; this allows only for fiber-preserving varia- 
tions, i.e., if 4>{X) S Y x and is a new section, then (j> 6 Y x . In general, one 
should allow for arbitrary variations in TY, when <p € Y x for some i/i, 



10 



Introducing a splitting of the tangent bundle into a vertical and a horizontal 
parts, T y Y = V y Y ©H y Y (H y Y is not uniquely defined), one can decompose 
a general variation into a vertical and horizontal components, respectively. 

Explicit calculations show (see Marsden, Pekarsky, Shkoller, and West 
[2000]) that while both vertical and arbitrary variations result in the same 
Euler-Lagrange equations, the Cartan form obtained from the vertical vari- 
ations only is missing one term (corresponding to the d n+l x from on X); 
the horizontal variations account precisely for this extra term and make the 
Cartan form complete. 

One can account for general variations either by introducing new "tilted 
sections" , or by introducing some true new sections that compensate for the 
horizontal variation. The later can be implemented in the following way. 
Let U C X be a smooth manifold with smooth closed boundary. Define the 
set of smooth maps 



For each p G C, set p>x = kxy V 7 an d Ux = t^xy ° f(U), so that (px '■ U — ► 
Ux is a diffeomorphism and <p o tp x is a section of Y. The tangent space 
to the manifold C at a point tp is the set T V C defined by 

{V G C°°{X, TY)\tty,ty oV = <pk Tttxy ° V = V x , a vector field onX}. 

Arbitrary (i.e., including both vertical and horizontal) variations of sec- 
tions of Y can be induced by a family of maps (f defined through the action 
of some Lie group. Let Q be a Lie group of ttxy bundle automorphisms rjy 
covering diffeomorphisms rjx- Define the action of Q on C by composition: 
rjY ■ ip = t]y o tp. Hence, while p o tp^ 1 is a section of t^u x X > Vy • f induces a 
section r] Y o (<p ■ p" 1 ) o t?" 1 of vr r;x(c/x)i y . 

A one parameter family of variations can be obtain in the following way. 
Let e i— > 7?y be an arbitrary smooth path in Q with r/y = e, and let V G T^C 
be given by 



Define the action function 



C = {tp:U 



Y\ttxy o tp :U ^> X\s an embedding}. 




and call p a critical point ( extremum ) of S if 

dS{p) ■ V = ^(r/y ■ yj) = 0. 



e=0 
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The Euler-Lagrange equations and the Cartan form can be obtained by 
analyzing this condition. We summarize the results in the following theorem 
from Marsden, Patrick, and Shkoller [1998] which illustrates the application 
of the variational principle to multisymplectic field theory. 

Theorem 2.1. Given a smooth Lagrangian density C : J^Y — > A n+1 (X) ; 
there exist a unique smooth section D EL C G C°°(y", A n+1 (X) ® T*Y)) 
(Y being the space of second jets of sections) and a unique differential form 
@c £ A n+1 (J 1 y) such that for any V £ T^C, and any open subset Ux such 
that U x n dX = 0, 

dS(ip) ■ V = [ D EL C(j 2 (<p o p" 1 )) ■ V + / j\<p o ^ 1 )*b' 1 (V) J Qc\. 

JUx JdU x 

(2.7) 

Furthermore, 

D EL C(j 2 (^o Vx 1 )).V=j 1 ^o i p x 1 y[j 1 (V)jn c ] inU x . (2.8) 

In coordinates, the action of the Euler-Lagrange derivative D E lC on Y" is 
given by 

r dL d 2 L 

d el cu\<p o vx 1 )) = [oj;UH<p o vx 1 )) - g^rU 1 ^ ° vx 1 )) 



dy b dv 



a 



d 2 L 

■ij X i}P 0( Px)) ' (<P 0i Px) 



dy a A d n+i x, (2.9) 



w/u/e the form Qc matches the definition of the Cartan form, obtained via 
Legendre transformation and has the coordinate expression 

BT f FIT \ 

&c = ^rdy" A d% + (L - _«- M J d"* 1 *. (2.10) 



Corollary 2.1. T/ie (n + \)-form 0£ defined by the variational principle 
satisfies the relationship 

for all holonomic sections 3 € C°°(ttx j^y)- 
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Another important general theorem, which we quote from Marsden, 
Patrick, and Shkoller [1998], is the so-called multisymplectic form for- 
mula 



Theorem 2.2. If cf) is a solution of the Euler- Lagrange equation fi2.9j ), then 
(/(^Pi))* [j 1 VJj 1 WjSl c ]=0 (2.11) 



'dUx 

for any V, W which solve the first variation equations of the Euler- Lagrange 



equations, i.e. any tangent vectors to the space of solutions of (2.£>). 



This result is the multisymplectic analog of the fact that the time t map 
of a mechanical system consists of canonical transformations. See Marsden, 
Patrick, and Shkoller [1998] for the proofs. 

Finally we remark that in order to obtain vertical variations we can 
require fx (and, hence, V 9 ^ 1 ) *° ^ e * ne identity map on X. Then, (f> = (poip^ 1 
becomes a true section of the bundle Y . 

Euler-Lagrange Equations. Treating (J 1 ! 7 , Qc) as a multisymplectic 
manifold, the Euler-Lagrange equations can be derived from the following 
condition on a section <f> of the bundle Y: 

(iV)*(>vjfi £ ) = o ) 

for any vector field W on J^Y (see Gotay, Isenberg Marsden [1997] for the 
proof). This translates to the following familiar expression in coordinates 

J^V)-^(^-V))=o, (, 12) 



which is equivalent to equation (2.9). 



Substituting the Lagrangian density fl2.3|) into equation ( |2.12| ) we obtain 
the following Euler-Lagrange equation for a continuous medium: 

d ( dW 




dW dg bc 4 



where 
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is the covariant time derivative, which corresponds to material accelera- 
tion, with 

c = 1 cd ( d 9ad d%d _ dg ab 
lab 2 9 V dy b dy a dy d 

being the Christoffel symbols associated with the 'field' metric g. We remark 
that all terms in this equation are functions of x and t and hence have the 
interpretation of material quantities. 

Equation ( |2.13 ) is a PDE to be solved for a section 4>(x,t) for a given 



type of potential energy W. As the gravity here is treated parametrically, 
the term on the right hand side of ( [2.13 ) can be thought of as a derivative 



with respect to a parameter, and we can define a multisymplectic analogue 
of the Cauchy stress tensor a as follows 

a ah = 2 4^{j 1 4>) : X^R, (2.14) 



J dg 



ab 



where J = det [F] \J det [g] / det [G] is the Jacobian. Equation ( |2.14| ) is known 
in the elasticity literature as the Doyle-Ericksen formula (recall that our p 
corresponds to pn e f, so that the Jacobian J in the denominator disappears). 
Another important remark is that the balance of moment of momentum 



T 

a = a 



follows from definition ( 2.14| ) and the symmetry of the metric tensor g. 



Finally, in the case of Euclidean manifolds with constant metrics g and 
G, equation ( |2.13| ) simplifies to 

d 2 (/) a d f dW ,. 1 



Barotropic Fluid 



"1# " a? ( W iV) 1 ■ (2 ' 15) 



For standard models of barotropic fluids, the potential energy of a fluid 
depends only on the Jacobian of the fluid's "deformation", so that W = 
W(J(g,G,v)). For a general inhomogeneous barotropic fluid, the material 
density is a given function p(x). In material representation, this formalism 
also includes the case of isentropic fluids in which there is a possible depen- 
dence on entropy. Since, in that case, entropy is advected, this dependency 
in the material representation is subsumed by the dependency of the stored 
energy function on the deformation gradient. 

2 In spatial representation, of course one has to introduce the entropy as an independent 
variable, but this naturally happens via reduction. See Holm, Marsden and Ratiu [1998] 
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The Legendre transformation can be thought of as defining the pressure 
function P. Notice that 



Pa 



8W r——- 

-/te\/det[G] 



dW dJ r——- 



dv a i dJ dv . 

and define the pressure function to be 

dW 



dW i i , 



P{cj>,x)=-p{x)-^j{j l 4>{x)) : CxX 



(2.16) 



Then for a given section <j), P(4>) : X — > R has the interpretation of the 
material pressure which is a function of the material density. In this case, 
the Cauchy stress tensor defined by ( 2.1 4| ) is proportional to the metric with 
the coefficient being minus the pressure itself: 



a ab (x) 



2pdW dJ , 2P 1 o6 4 

^"^T« — V 9) = — r J o5 U ( 
J dJ dg ab J 2 



-P(x) g ab (y(x)). 



We remark that this can be a defining equation for the pressure from which 
( |2.16| ) would follow. With this notation the left hand side of the Euler- 
Lagrange equations ( 2.13 ) becomes 



P9ab 



P9ab 




+ (I) + {II) 



+ 2 J9 o r ' 



(2.17) 



where terms (I) and (II) arise from differentiating det[v] and (v 1 ) a k and 
cancel each other. The right hand side of ( 2.13 ) is given by 



-P 



dW dg bc 
dgbc dy a 



-p 



dW dJ dg 



be 



dJ dg bc dy a 



P T _.bc 9 9bc 

2 J9 dy"' 



Notice that the last term in ( |2.17| ) and in the equation above coincide, so 
that the Euler-Lagrange equations for the barotropic fluid have the following 

for related results from the point of view of the Euler-Poincare theory with advected 
quantities. 
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form 

where the pressure depends on the section (p and the density p and is de- 



fined by ( 2.16 ). Both the metric g a b and the Christoffel symbols 7^ fe in the 
covariant derivative are evaluated at y = (p(x,t). 

One can re- write ( |2.18| ) introducing the spatial density p sp = p/J 
and defining the spatial pressure p(y) by the relation P(x) = p(y(x)) = 
p(4>t(x). This yields 

DaV , . 1 , . , 

-f-r(x,t) = gradpo (x,t), 

Dt p sp 

where V = <j). Compare this to the equations for incompressible ideal hy- 
drodynamics in §[|. 

Elasticity 

The Legendre transformation defines the first Piola-Kirchhoff stress tensor 
Va- It is given, up to the multiple of — 1/ y det [G] , by the matrix of the 
partial derivatives of the Lagrangian with respect to the deformation gradi- 
ent: 

dW 1 

PSfax) = pW^tfftx)), (2.19) 

and for a given section Va is a stress tensor defined on X. 

Notice that the first Piola-Kirchhoff stress tensor is proportional to the 
spatial momenta, V a % = —p a % I sj det [G] . The coefficient yj det [G] arises from 
the difference in the volume forms used in standard and multisymplectic 
elasticity. In the former, the Lagrangian density is integrated over a space 
area using the volume form pq = w det[G]d n x associated with the metric 
G, while in the latter, the integration is done over the space-time using 
d n+1 x = dt A d n x. We also remark that though traditionally the first Piola- 
Kirchhoff stress tensor is normally taken with both indices up, our choice 
is more natural in the sense that it arise from the Lagrange transformation 
( |2.19| ) which relates TV with the spatial momenta. 

Using definitions ( 2.14| ) and fl2.19|) , we can re- write equation ( 2.13| ) in 



the following form 

P9ah (^4) = V a \ + 7 2c(*V> c i " J9M° dc ), (2.20) 
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where we have introduced a covariant divergence according to 



V a \ = BWV = + VjV% - V^lF^. 

Here T % - k are the Christoffel symbols corresponding to the base metric G on 
B C X (see, for example, Marsden and Hughes [1993] for an exposition on 
covariant derivatives of two-point tensors). 

We emphasize that in ( |2.20| ) there is no a-priori relationship between the 
first Piola-Kirchhoff stress tensor and the Cauchy stress tensor, that is, W 
has the most general form W(x,G, g, v). Such a relationship can, however, 
be derived from the fact that for standard models of elasticity the stored 
energy function W depends on the deformation gradient F (i.e. on v) and 
on the field metric g only via the Green deformation tensor C given by ( |2.2| ) , 
that is W = W(C(v, g)). Thus, the partial derivatives of W with respect to 
g and v are related, and the following equation 

TV = J(<?F- l ) a l 



follows from definitions ( 2. 14| ) and (|2.19| ). This relation immediately follows 



from the form of the stored energy function; it recovers the Piola transfor- 
mation law, which in conventional elasticity relates the first Piola-Kirchhoff 
stress tensor and the Cauchy stress tensor. Substituting this relation in 
Q2.20p one easily notices that the last term on the right hand side cancels, 
so that the Euler-Lagrange equation for the standard elasticity model can 
be written in the following covariant form 

p ^- = mvv, (2.21) 

where V = <j>. For elasticity in a Euclidean space, this equation simplifies 
and takes a well-known form: 

dV a dV ai 
P- 



dt dx % 



3 Constrained Multisymplectic Field Theories 

Multisymplectic field theory is a formalism for the construction of Lagran- 
gian field theories. This is to be contrasted with the formalism in which 
one takes the view of infinite dimensional manifolds of fields as configura- 
tion spaces. The multisymplectic view makes explicit use of the fact that 
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many Lagrangian field theories are local theories, that is, the Lagrangian 
depends only pointwise on the values of the fields and their derivatives. In 
formulating a constrained multisymplectic theory, we will therefore only be 
concerned with the imposition of pointwise constraints $(7), 7 £ J 1 !', de- 
pending on point values of the fields and their derivatives. In the current 
work we also restrict our attention to first-order theories, in which only first 
derivatives of the fields are considered. 

Despite the pointwise nature of the Lagrangian £(7), 7 G J X Y ■, the 
variational principle assumes variations of local sections over some region 
U C X, that is, it is the action S((j)) = J' [/ £(j 1 0) as a function of sections 
that is being minimized. In order to use the theory of Lagrange multipliers 
to impose the constraints, it is therefore necessary to form a function ^f((j>) of 
local sections which is defined through point values of the constraint ^(j 1 ^) 
evaluated at the first jets of sections. It is then possible, however, to use 
the pointwise nature of the Lagrangian and the constraint function to derive 
a purely local condition, the Euler-Lagrange equations, for the constrained 



field variables. We will make these ideas precise in Section 3.2 



For holonomic constraints it is well known that Hamilton's principle 
constrained to the space of allowable configurations gives the correct equa- 
tions of motion. Hamilton's principle can be naturally extended by either 
extremizing over the space of motions satisfying the constraints (so-called 
vakonomic mechanics) , which is appropriate for optimal control, but not for 
dynamics, or by requiring stationarity of the action with respect to varia- 
tions which satisfy the constraints (the Lagrange-d'Alembert or virtual work 
principle). The equations of motion derived in each case are, however, dif- 
ferent. 

Derivations from balance laws (Jalnapurkar [1994]), evidence from ex- 
periments (Lewis and Murray [1995]) and comparison to Gauss' Principle of 
Least Constraint and the Gibbs-Appell equations (Lewis [1996]) indicates 
that it is the Lagrange-d'Alembert principle which gives the correct equa- 
tions of motion; see Bloch and Crouch [1999] for further discussion and 
references. 

While the subject of linear and affine non-holonomic constraints is rel- 
atively well-understood (see Bloch, Krishnaprasad, Marsden, and Murray 
[1996]), it is less clear how to proceed for non-linear non-holonomic con- 
straints. Part of the problem lies in the lack of examples for which the 
correct equations are clear from physical grounds. In the context of con- 
strained field theories, however, there are many cases where nonlinear con- 
straints involving spatial derivatives of the fields need to be applied, such 
as incompressibility in fluid mechanics, and it is clear what the physically 
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correct equations should be. Here we deliberately avoid the use of the term 
non-holonomic, to avoid confusion with its standard meaning in the ODE 
context, where it applies only to time derivatives. Other examples of nonlin- 
early constrained field theories include constrained director models of elastic 
rods and shells. 

The fact that the constraints involve only spatial and not time deriva- 
tives means that imposing the constraints is equivalent to restricting the 
infinite-dimensional configuration manifold used to formulate the theory as 
a traditional Hamiltonian or Lagrangian field theory. In this case, the con- 
straint is simply a holonomic or configuration constraint and it is known 
that restricting Hamilton's principle to the constraint submanifold gives the 
correct equations for the system. 

3.1 Lagrange Multipliers 

The Lagrange multiplier theorem naturally makes use of the dual of the 
space of constraints. In a finite-dimensional setting this is a well defined 
object, with all definitions being equivalent. When considering infinite- 
dimensional constraint spaces, however, the issue of what is being used as 
the dual becomes less clear and more important. 

We shall consider constrained multisymplectic field theories for which 
the constraint space is the space of smooth sections of a particular vector 
bundle. In the case of the incompressibility constraint, the vector space 
is one-dimensional and the constraint bundle is, effectively, the space of 
real valued functions on the base space X. A dual of the constraint space 
is then defined with respect to an inner product structure on the vector 
bundle. This is made explicit in the following statement of the Lagrange 
multiplier theorem where we assume that fields and Lagrange multipliers 
are sufficiently regular (see Luenberger [1969]). 

Theorem 3.1 (Lagrange multiplier theorem). Let ttm,£ '■ £ ~ > M. be 
an inner product bundle over a smooth manifold Ai, Vl/ a smooth section 
of ttm,£; an d h : Ai — > R a smooth function. Setting Af = X I /_1 (0), the 
following are equivalent: 

1. (p £ Af is an extremum of h\j\f 

2. there exists an extremum (p £ 6 of h : 6 — ► K such that itme^) = <P 
where h(cp) = h(irM,e(<p)) ~ ( l P,^(^M,e( < P)))£- 
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If £ is a trivial bundle over A4, then in coordinates of the trivialization 
we have <p = (ip, A), where A:.M^£/.Misa Lagrange multiplier function. 

In the next section we shall use this theorem to relate the constrained 
Hamilton's principle with the extremum of the augmented action integral 
which contains the constraint paired with a Lagrange multiplier. Both of 
them result in constrained Euler-Lagrange equations. We shall furthermore 
demonstrate that, using the trivialization coordinates, these equations can 
be equivalently obtained from a Lagrangian defined on an extended con- 
figuration bundle. In this picture, the Lagrange multiplier corresponds to 
a new field, which extends the dimension of the fiber space, and the aug- 
mented Lagrangian contains an additional part corresponding to the pairing 
of this field with the constraint. The Euler-Lagrange equations of motion 
then follow from unconstrained Hamilton's principle in a standard way. 

3.2 Multisymplectic Field Theories 

In the setting above, the configuration bundle is a fiber bundle ttx y ■ Y — > X 
and 7Ty jxy '■ J l Y — > Y is the corresponding first jet bundle with and y a 
being a local coordinate system on X and Y respectively, and v a ^ the fiber 
coordinates on J 1 Y '. 

Choose the configuration manifold Ai to be the space C of smooth sec- 
tions cf> of ttx,y- R eca h that for a Lagrangian density C : J l Y — > A n+1 X, a 
section <j) E Ai is said to satisfy Hamilton's principle if <fi is an extremum of 
the action function S((p) = J x C(j l 4>) : Ai — > M. Choose the h above to be 
the action function S and use S instead of h. 

To apply the Lagrange multiplier theorem we need to define constraints 
as a section of some bundle £ — > A4 (below called the constraint bundle). 
As mentioned above, we restrict our attention to constraints $ which de- 
pend only on point values of the fields and their derivatives. Using such 
constraints we can construct induced constraints ^ according to (|3.1| ). This 
is made precise below. We point out, however, that our treatment excludes 
inherently global constraints, such as those on the inverse Laplacian of the 
field, which can not be derived from pointwise values. 

On the other hand, we also exclude from the consideration a (simple) 
case when the constrained subbundle of J l Y can be trivially realized as the 
first jet of some subbundle of Y. 

Define an inner product vector bundle ttx,v ■ V — > X with the inner 
product denoted by (■, -}y whose fibers are isomorphic to M. n . Let C°°(V) be 
the inner product space of smooth sections of ttx,v with the inner product 
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given by 

(a, b) = [ (a(x),b(x)) v d n+l x. 
Jx 

The constraint function is an M. n - valued function on J^Y: 

$ : J l Y -» R n . 

We say that a point 7 € J^Y satisfies the constraint if $(7) = 0. By 
restricting 3> to the space of first jets of sections (ft of Y, we can define the 
induced constraint function ^ from $ by setting 

*(0)(aO = (3.1) 

for all € and x € X. By construction, ^ is a map from the space A4 of 
sections of Y to the space C°°(V) of sections of V, hence it can be thought of 
as a smooth section ^ : A4 — > £ of the constraint bundle £. This bundle 
is the trivial inner product bundle given by M x C°°(V) over j\A. Then, a 
configuration (ft € .M is said to satisfy the constraints if •3>((j 1 ^)(x)) = for 
all x £ X, that is, the section *$>((ft) must be a zero function on X. This 
implies that the space of configurations which satisfy the constraints is given 
by j\f = V- 1 (0). 

The constrained Hamilton's principle now seeks a (ft £ J\f which is an 
extremum of S\j\f. The Lagrange multiplier theorem given in the previous 
subsection can be applied to conclude that this is equivalent to the existence 
of (ft € £ with ttm,£(4>) = 4> which is an extremum of S. Using the coordinates 
of the trivialization of £ we can write <ft = ((ft, A), where (ft = 7Tm,s(^) is the 
base point, i.e. section (ft of Y, and A is thought of as a section of the bundle 
ttx,v, i-e. an K n -valued function on X. Then S : £ — > R is given by 

S$) = S(cft)-(\,*(cft))e 

= [ L{(j l cft){x))cr +1 x- ( (A(x),cD(( J V)(x)))v^ +1 x 
Jx Jx 

= [ [L(C/V)(z)) " (\(x),m l <P)(x)))v]d n+1 x 
Jx 

In the next section, we demonstrate these constructions for the incom- 
pressibility constraint for continuum theories. 

The requirement that S be stationary with respect to variations in A at 
the point (ft implies that 
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for all variations <5A, and thus that <I > ((j 1 </>)(x)) = for all x € X. This 
therefore recovers the condition that 4> must satisfy the constraints. 

Stationarity of S with respect to variations in <fi can be used to derive 
the constrained Euler- Lagrange equations, which have the form 

£(£<«V)M))-J>-'«M) 

+ J^(O'V)M)) - £: (aw, |!r(« V) (*))) = °- M 

Alternatively, one can handle the constraints by introducing another 
bundle, denoted by E, which is a product bundle over X with fibers diffeo- 
morphic to Y x x V x - One can think of E as a configuration bundle of the 
corresponding unconstrained system whose Lagrangian contains an addi- 
tional part corresponding to the pairing of the constraint with the Lagrange 
multiplier: 

L$ = L + (A,$) v . 

The Euler-Lagrange equations of motion then follow from unconstrained 
Hamilton's principle in a standard way and coincide with ( |3.2[ ). We work 
out the details for the incompressibility constraint in the next section. 

4 Incompressible Continuum Mechanics 

In this section we shall consider the incompressibility constraint using the 
multisymplectic description of continuum mechanics. The main issue is a 
proper interpretation of the constraint using the Lagrange multiplier for- 
malism developed in the previous section. 

4.1 Configuration and Phase Spaces 

Here we briefly summarize the results. See the analogous parts of §|2| for 
more details. 

Extended Covariant Configuration Bundle. The fibers of V in this 
case are one-dimensional and sections eft = (<f>, A) of E contain both the 
deformation field and the Lagrange multiplier, i.e., E denotes a bundle over 
X whose fibers are diffeomorphic to the product manifold Mxl with the 
projection map 

ttxe ■ E -> X; (x,t,y, A) h-> (x,t). 
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Here A is a section of the trivial bundle X x K over X, which can be thought 
of as a function \(x, t) on X. The phase space is then the first jet bundle J X E 
with coordinates 7 = (x^, y a , A, u a M , (3^); the first jet extension of a section 
4> = (4>, A) has the following coordinate representation (x' 1 ,y a , A, d^cj) 11 , d^X). 

The Dual Jet Bundle. We can consider the affine dual bundle J l E* 
as a "vertically invariant" subbundle Z of the bundle A = A n+1 E of all 
(n + l)-forms on E. Elements of Z can be written uniquely as 

z = Ud n+1 x + p a ^dy a A d n Xf , + TT^dX A d n Xfl , 

where d n x }1 = d^d n+1 x, so that (x^,y a , A, II, p a ^, tt^) give coordinates on Z. 

The canonical (n + l)-form is constructed in a standard manner and in 
the above coordinates has the following representation 

6 = p a "dy a A d n x^ + vr^A A d n Xfl + Ud n+1 x. 

We set n = -d9. 

The primary constraint manifold £ is a subbundle of the dual jet bundle 
and corresponds to the incompressibility constraint. The pull-back of the 
inclusion map i<r : £ — > J 1 ^* defines a degenerate (n + 2)-form on £. We 
shall discuss the explicit form of the constraint in the next subsection. 



Incompressibility Constraint. Recall that such a constraint in, for ex- 
ample, incompressible fluid dynamics, is a reflection of the divergence-free 
property of the Eulerian fluid velocity and, hence, has a pointwise character. 
The divergence-free character of the velocity field arises from the require- 
ment that the particle placement map be volume preserving at each instant 
of time. Then, according to the general theory of constrained multisymplec- 
tic fields outlined above, it can be obtained from a pointwise constraint <3? 
defined on the first jet bundle J l Y . 

For 7 = (x tM ,y a ,v a ^ t ) £ J l Y we impose the constraint $(7) = on the 
Jacobian of the deformation, where 



where we have used the definition of J given in (2.1). Restricting $ to the 



first jet of a section cj) results in a constraint on the matrix of spatial partial 



derivatives djc 
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For the Lagrange multiplier itself, we choose the following Ansatz 

A(x) = ^/det[G]P(x) : X -> M, (4.2) 

where P will be shown later to have the interpretation of the material pres- 
sure. Equation (|4.2j ) guarantees that A transforms like a density under the 
transformations of the base manifold X, so that the pairing of A and <£, 
defined by integrating over X, has the correct transformation law. 



4.2 Lagrangian Dynamics. 

As we have already mentioned, the main distinguishing feature of incom- 
pressible models of continuum mechanics is the presence of the constraint 
(4.1). We shall now explain how this modification to the Lagrangian alters 
the Legendre transform as well as the Euler-Lagrange equations. 



The Lagrangian Density. The Lagrangian density C : J E —* A n+ X 
for the multisymplectic model of incompressible continuum mechanics is a 
smooth bundle map defined by 

£$(7) = (L( 7 ) + A • $(7)) d n+1 x = K - P + A • $ d n+l x, (4.3) 

where L (i.e. K and IP) is given by ( [2.3] ) and depends on the choice of the 
stored energy function W. 



The Legendre Transformation. For the above choice of the Lagrangian, 
the Legendre transform thought of as a fiber preserving bundle map FC$ : 
J l E — > J l E* over E is degenerate due to the constrained character of 
the dynamics. Indeed, the Lagrange multiplier A is a cyclic variable as 
the Lagrangian ( |4.3j ) does not depend on its derivatives, [3^. Hence, the 
conjugate momentum to A is identically zero: 7r M = dL$/df3^ = 0. The set 
{-7r M = 0} defines the primary constraint set as a subset of the dual bundle 
J 1 £7* to which we restrict the Legendre transformation to make it non- 
degenerate. The rest of the momenta are given by the following expressions 

n <9L$ 



Pa 



dv a o 



P9abV 



v/det[G], 



Pa 



(4.4) 



n 



1 



a „,b 



dW 



pW - P ( J{n - 1) + 1) 



v/detM- 
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Euler-Lagrange Equations Using the trivialization (0, A), we now con- 
sider the Euler-Lagrange equations for a section (ft of E, both with respect 
to the deformation <p and with respect to the Lagrange multiplier A. The 
former can be written in coordinates as follows 

|^W)=0. (4.5) 
dy a u YJ dx» \dv a ^ J YJ J y ' 

The Euler-Lagrange equation for A trivially recovers the constraint $ = 
itself restricted to the first jet : 

^ (i 1 ^) " 7^7 ( ^U^)) = Hj^W+'x = (J(iV) - l)d n+1 x = 0. 



(4.6) 

These two equations are to be solved for the Lagrange multiplier A (equiv- 
alently, for the pressure P) and for the section </>. 

Substituting Lagrangian ( |4.3| ) into ( |4.5| ), we obtain the Euler-Lagrange 
equation ( p. 13 ) modified by the pressure term: 



v ^t\G]dx k 

dW dg bc , &P fc ! 

Notice that in the case of parameterized non-constant metrics, the extra 
pressure term in ( |4.4| ) gives rise to the term 



which follows from the chain rule applied to d x j g(y(x)). This term exactly 
cancels another term coming from differentiating the constraint with respect 
to y due again to the composition g = g{y) : 

a|^ = ^jg^v^m 

dy a 2 dy a v L J 

and other cancellations occur as in equation ( 2,17| ), 

In the case of Euclidean manifolds with constant metrics g and G, the 
Euler-Lagrange equations simplify to 

d 2 6 n 8 ( dW , .i .A dP 



^ dt 2 dx k \dv a k yJ ' r ' ) dx k 
together with the constraint (|L6|). 



y i \ or i k 1 

- ^(O.* ^V) (4-8) 
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4.3 Incompressible Ideal Hydrodynamics 

For fluid dynamics, the stored energy term in the Lagrangian is a constant 
function precisely because of the incompressibility constraint. Indeed, as we 
have mentioned above, W in ideal fluid models is a function of the Jacobian 
J, but the latter is constrained to be 1. For simplicity, consider an ideal 
homogeneous incompressible fluid, so that the material density p is constant, 
and we can set p = 1 (for inhomogeneous fluids the dependence of material 
density on the point x is implicit in the pressure function P). 

The Lagrangian is given by (|4.3|) with F = const. Hence, two terms in 



equation (4.7) which correspond to the derivatives of W vanish, so that the 



dynamics of an incompressible ideal fluid is described by 




together with the constraint 



where we have used the fact that g = G. 

Compare (|4.9| ) with fl2.18|) and notice that the incompressibility con- 
straint J{j l (f)) = 1 implies that the spatial density p sp = p/J is constant, 
e.g., 1. Introducing the spatial pressure p = P o , the above equation 
can be written as 



Dt 



(x,t) = — gradp o cp(x,t), (4-11) 



where we have set p sp = 1. We remark again that the covariant derivative 
is evaluated at y = tp(x,t). 

A New Look at the Pressure. Here we shall demonstrate that the same 
equations of motion are obtained if the potential energy in the Lagrangian 
Q4.3|) is not set to a constant, but rather is treated as a function of the 
Jacobian, W = W(J). This will also clarify the relation between the two 
definitions of pressure that we have thus far examined. 

Recall the definition of the pressure function for barotropic fluids given 
by ( 2.16| ) as a partial derivative of the stored energy function W with respect 



to the Jacobian J. Compare this to the definition Q4.2j) of the pressure as a 
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Lagrange multiplier corresponding to the incompressibility constraint (4.1) 
(modulo a \J det [G] term). In this subsection we shall denote these objects 
by Pw and P\, respectively: 

Pw = -P^tti -Pa 



dJ' VdeTlGj 

The resulting Euler-Lagrange equations can be obtained by combining 
3. 18f) with ( |4.9| ) and are given by: 



Sab 



'zV 6 V_ dQFV + P A ) 




Dt J dx k 
together with the constraint fl4.10| ). We can define a new pressure function 

P = P W + P\- (4.12) 

Notice that when the constraint J = 1 is enforced by the Euler-Lagrange 
equation ( [4.6D , P\y(J) = const, so that P = P\ + const. This is equivalent 
to a re-definition of the Lagrange multiplier A. At the same time, the above 
Euler-Lagrange equation coincides with (|4.9| ) because d^P = dk^x- 



Relation to Standard Ideal Hydrodynamics. Recall the Lie-Poisson 
description of fluid dynamics as a right-invariant system on the group T>^(M) 
of volume-preserving diffeomorphisms of a Riemannian manifold (M, G) . 
Here we follow Marsden and Ratiu [1999] and Arnold and Khesin [1998], 
using our notations. The Lie algebra of V^(M) is the algebra of divergence- 
free vector fields on M tangential to the boundary with minus the Jacobi-Lie 
bracket. The L 2 inner-product on this algebra is given by 

(u,v) L 2 = / (u(x),v(x)) G n, 
JM 

where is the Riemannian volume form on M. 

We extend this metric by right invariance to the entire group. The 
resulting Riemannian manifold with right invariant 1? metric, denoted by 
(2}^(M), L 2 ), is the configuration space for the Lie-Poisson or Euler-Poincare 
model of ideal hydrodynamics. Its tangent bundle is the phase space, so that 
(rjt,fjt) are the basic "coordinates"; here rjt £ T>^(M) is a diffeomorphism 
that transforms the reference fluid configuration to its configuration at time 
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t. Then, using the kinetic energy of fluid particles as a Lagrangian, one 
obtains the following covariant equations of motion: 



Df) 
~Dt 



(x) = — gradp o 7](x), 



(4.13) 



where 



Dfj 
~Dt 



denotes covariant material time derivative with respect to the metric (T v 
denotes the connection associated to the metric) and p is the spatial pressure. 
Notice that covariant derivative is evaluated at rj(x). 

Now define rjt(x) = rj(t,x) to be the flow of the time-dependent vector 
field u(t,x), so that dtrj(t,x) = u(t, r)(t, x)) . Then composing ( 4.13 ) on the 
right with r/ _1 gives the classical Eulerian description of incompressible ideal 
fluids : 

dtu(t, x) + (u ■ V)u = — gradp, div u = 0. 

Taking the divergence of both sides of this expression yields the equation 
for the pressure 



Ap = - div ((it • V)u) 



(4.14) 



One readily notices that equations (4.11) and ( 4.1 3| ) coincide provided 
r]t{x) = <fi(x,t). Upon this identification, the Euler-Lagrange equations for 
the multisymplectic model of incompressible ideal hydrodynamics recover 
the well-known evolution of fluid diffeomorphisms ( 4.13| ). Similarly, taking 
the divergence of both sides of ( 4.1 1) ) results in the Poisson equation on the 
pressure ( |4. 14 ) . 



4.4 Incompressible Elasticity 

In a manner similar to the previous subsection, we modify the elasticity 
Lagrangian by the constraint and extend the phase space to include the 
Lagrange multiplier. Recall that the stored energy is a function of the Green 
deformation tensor W = W(C) and use the definition of the first Piola- 
Kirchhoff stress tensor V a l (|2.1S| ) to write down the equations of motion: 
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together with the constraint ( |4.6| ). The above equation can be written in a 
fully covariant form 



where V = 4> is the velocity vector field, V is the first Piola-Kirchhoff stress 
tensor, and p is the spatial pressure. 



5 Symmetries, Momentum Maps and Noether's 
Theorem 

We already mentioned in the introduction that homogeneous fluid dynamics 
has a huge symmetry, namely the particle relabeling symmetry, while stan- 
dard elasticity (usually assumed to be inhomogeneous) has much smaller 
symmetry groups, such as rotations and translations in the Euclidean case. 
While inhomogeneous fluids (especially the compressible ones) are of great 
interest, the results worked out in §|0] only apply to homogeneous fluid 
dynamics, when the symmetry group is the full group of volume-preserving 
diffeomorphisms P«. However, these results can be generalized to inhomo- 
geneous fluids, in which case the symmetry group is a subgroup C T>^ 
that preserves the level sets of the material density for barotropic fluids, or 
a subgroup Df l ' ent c T>^ that preserves the level sets of the material density 
and entropy for isentropic fluids. 

A general model of continuum mechanics will have the metric g isometry 
as its symmetry. In particular, the group of rotations and translations is a 
symmetry for models of fluid dynamics and elasticity in Euclidean spaces. 
The later is treated in Marsden, Pekarsky, Shkoller, and West [2000], where 
the overall emphasis is on continuum mechanics in Euclidean spaces. 

The only symmetry which is universal for non-relativistic continuum 
mechanics is the time translation invariance. This is due to the fact that the 
base manifold is a tensor product of the spatial part and the time direction, 
rather than a space-time, so that all material quantities, such as density 
p, metric G, etc. depend only on x € B C X. In this section we shall 
treat these symmetries separately. We start with the particle relabeling 
symmetry, introducing the necessary notations. 
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5.1 Relabeling Symmetry of Ideal Homogeneous Hydrody- 
namics 

In this subsection we shall consider both the barotropic model and the in- 
compressible model of ideal homogeneous fluids with fixed boundaries at the 
same time. Their corresponding Lagrangians differ only by the constraint 
term and both are equivariant with respect to the action of the group of 
volume-preserving diffeomorphisms. 

The Group Action. The action of the diffeomorphism group T>^{B) on 
the (spatial part of the) base manifold B C X captures precisely the meaning 
of particle relabeling. For any r/ € V^(B), denote this action by rjx '■ (x, t) i— > 
(r](x),t). The lifts of this action to the bundles Y and E are given by 
rjy : (x,t,y) i-> (rj(x),t,y) and rj E : (x,t,y,X) h-> (j](x),t, y, A), respectively. 
Both lifts are fiber-preserving and act on the fibers themselves by the identity 
transformation. The coordinate expressions have the following form: 

rf 3 x =ld-t,rf x =rf{x), r, Y = 5 a b y\ r,% = («5 6 V, Id- A). (5.1) 

Jet Prolongations. The jet prolongations are natural lifts of automor- 
phisms of Y to automorphisms of its first jet J l Y and can be viewed as 
covariant analogues of the tangent maps (see Gotay, Isenberg and Marsden 
[1997]). 

Let 7 be an element of J^Y and 7 be a corresponding element of the 
extended phase phase space in coordinates 7 = (x M , y a , v°) and 7 = 

(x^, y a , A, v^, Pn). The prolongation of rjy is defined by 

Vj^y (7) = Ti] Y 07 oT^ 1 , t)jie(i) = T VE ojoTrjx 1 . (5.2) 

We shall henceforth consider r/jig, since it includes rjjiy as a special case. 
In coordinates, we have 

W7) = (W), t; y\ A; v\, v a m ((f^)" 1 ) " ',Po,Pm ■ 

If £ is a vector field on E whose flow is rj e , then its prolongation j 1 ^ is the 
vector field on J 1 !? whose flow is j 1 ^), that is j 1 ^ oj 1 ^) = (d/de)j 1 (r/ e ). 
In particular, the vector field £ corresponding to t]e given by ( |5.1|) has 
coordinates (£ l , 0,0,0) and is divergence-free; its prolongation j 1 ^, which 
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corresponds to the prolongation rjji E of 1]e, has the following coordinate 
expression: 

i^e,0;0,0;0,-v a m ^-;0,-f3 m ^-j . (5.3) 

Noether's Theorem. Suppose the Lie group Q acts on C and leaves the 
action S invariant. This is equivalent to the Lagrangian density C being 
equivariant with respect to Q, that is, for all r\ £ Q and 7 £ ^Y, 

£(^y(7)) = (Vx 

where (n^ 1 )* £(7) = {tix)*£>{h) is a push-forward; this equality means equal- 
ity of (n + l)-forms at r)(x). Denote the covariant momentum map on 
J X Y by Jc £ L(q, A n (J 1 Y)). It is defined by the following expression 

j 1 (0^n jC = dJ c (0 (5.4) 

and can be thought of as a Lie algebra valued ra-form on ^Y. 

Recall that is a solution of the Euler-Lagrange equations if and only if 

(iV)*(WJfir)=0 

for any vector field W on J l Y . In particular, setting W = j 1 (£) and applying 
(j 1 ^)* to (|5.4|), we obtain the following basic Noether conservation law: 

Theorem 5.1. Assume that group Q acts on Y by irxY-bundle automor- 
phisms and that the Lagrangian density C is equivariant with respect to this 
action for any 7 £ J^Y . Then, for each £ £ q 

d((jV)*4c(0) =0 (5-5) 

/or any section <j> ofirxY satisfying the Euler-Lagrange equations. The quan- 
tity {j l 4>)* Jc{Q is called the Noether current. 

See Gotay, Isenberg and Marsden [1997] for a proof. 



The Variational Route to Noether's Theorem. The variational route 
to the covariant Noether's theorem on J l Y was first presented on pages 374- 
375 in Marsden, Patrick, and Shkoller [1998]. We shall briefly describe this 
formulation now. 

Recall the notations of the maps ip : U — ► Y and the corresponding 
induced local sections (p o ip~^ of Y from Subsection Here again it is 
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important to allow for both vertical and horizontal variations of the sections. 
Vertical variations alone capture only fiber preserving symmetries (i.e., spa- 
tial symmetries), while taking arbitrary variations allows for both material 
and spatial symmetries to be included. 

The invariance of the action S = J Ux £ under the Lie group action is 
formally represented by the following expression: 



S(t]y • <p) = S(tp) for all r/y £ Q. 



(5.6) 



Equation (|3j|) implies that for each r]y € Q, Tfy ■ (p is a solution of the 
Euler-Lagrange equations, whenever tp is a solution. We restrict the action 
of Q to the space of solutions, and let £e be the corresponding infinitesimal 
generator on C restricted to the space of solutions; then 



= (&JdS)fa) 



dU x 



fx 



b a (e)je £ ] 



Ux 



since the Lie derivative £jire\@c = by (5^) and Corollary 2.1 



Using (5.4), we find that f Ux d[j l {ip o <p~^ )*Jc(Q] = 0) anci since this 
holds for arbitrary regions Ux, the integrand must also vanish. Recall that 
(f> = ip o (p^ 1 is a true section of the bundle Y, so that this is precisely a 



restatement of the Noether's Theorem 5.1. 



Covariant Canonical Transformations. The computations of the mo- 
mentum map from definition (^J) can be simplified significantly in some spe- 
cial cases which we discuss here. A 7Txjiy-bundle map rjjiy '■ J l Y — ► J l Y 
covering the diffeomorphism rjx '■ X — > X is called a covariant canonical 
transformation ii rj jiyflc = It is called a special covariant canon- 
ical transformation if r]j lY @c = ©£• Recall that forms Qc an d @c can 
be obtained either by variational arguments or by pulling back canonical 
forms SI and from the dual bundle using the Legendre transformation WC 
From Gotay, Isenberg and Marsden [1997], any r\j\y which is obtained 
by lifting some action r\y on Y to J l Y , is automatically a special canonical 
transformation. In this case the momentum mapping is given by 

M0=j l ^®c- (5.7) 

We remark that the validity of this expression does not rely on the way 
in which the Cartan form was derived, i.e., for simplicity of the computa- 
tions in concrete examples, one can forgo the issues of vertical vs. arbitrary 
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variations in the variational derivation and obtain the Cartan form directly 
from the dual bundle by means of Legendre transformations. Then, evalu- 
ating this form on the prolongation of a vector of an infinitesimal generator 
gives the momentum n-form. 



Equivariance of the Lagrangian. To apply Theorem 5.1 to our case we 
need to establish equivariance of the fluid Lagrangians: 

Proposition 5.1. The Lagrangian of an ideal homogeneous barotropic fluid 



(2.5) and the Lagrangian of an ideal homogeneous incompressible fluid (4-5) 
are equivariant with respect to the T> fM (B) action ( \5. \ ): 

= te 1 )*^), 

for all 7 G J l E. 

Proof. First observe that the material density of an ideal homogeneous (com- 
pressible or incompressible) fluid is constant. Notice also that Lagrangians 
(2.3) and (4.3) differ only in the potential energy terms. Both these terms 



are functions of the Jacobian, which is equivariant with respect to the action 
of volume preserving diffeomorphisms given by (|5.3|). Indeed, 



/ (J (WTO)-«. = /^^g*t(„d*(g) j<r«* 

= (ixT (/(J(7)M" +1 x) 

due to the fact that det dirf = 1 for a volume preserving diffeomorphism rj; 
here / can be any function, e.g. the stored energy W or the constraint 



For the same reason, and the fact that (5.3) acts trivially on v a o, the 



kinetic part of both Lagrangians is also equivariant. □ 
Proposition 5.1 enables us to use (5.7) for explicit computations of the 



momentum maps for the relabeling symmetry of homogeneous hydrodynam- 
ics. We shall consider barotropic and incompressible ideal fluids separately 
because their Lagrangians and, hence, their momentum mappings are dif- 
ferent. 



Barotropic Fluid 

Using (fT?]) we can compute the Noether current corresponding to the rela- 
beling symmetry of the Lagrangian ( |2.3[ ) to be 
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fttTMO = QpW'V 6 -pw- pj\ y /dri\G\s le <rx k 



{9ab<P b <P%) pVtet\G\S k <rxo, (5. 
where j 1 ^ is the prolongation of the vector field £ and is given by (5.3). 



The differential of this quantity restricted to the solutions of the Euler- 



Lagrange equation is identically zero according to Theorem 5.1. Conversely, 
requiring the differential of (|5.8| ) to be zero for arbitrary sections <f> recovers 
the Euler-Lagrange equation. Indeed, computing the exterior derivative and 
taking into account that the vector field £ is divergence free, we obtain: 





which coincides with the Euler-Lagrange equation ( 2.1§| ). 



Incompressible Ideal Fluid 

Similar computations using Lagrangian ( [4.3[ ) with the potential energy set 
to a constant give the following expression for the Noether current corre- 
sponding to the relabeling symmetry: 

jHWMO = {^P9ab^ b ~ P^j Vd^G]ed n x k - 

{g ab ^ b ^ k ) pVdit[tffe*d»xo. (5.9) 



The assumptions of Theorem |5.1| are satisfied; hence the exterior dif- 
ferential of this Noether current d (j l (</))* Jcti 1 ^)) is equal to zero for all 
section (p which are solutions of the Euler-Lagrange equations. 

Now consider the inverse statement. That is, let us analyze whether 
the Noether conservation law implies the Euler-Lagrange equations for in- 
compressible ideal fluids. Computing the exterior differential of ( |5.9[ ) for an 
arbitrary section </> = {(f), A) we obtain: 





Here, we have used the fact that £ is a divergence free vector field on X. 
This is precisely the Euler-Lagrange equation (|4.9|) with the constraint J = 1 
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substituted in it. We point out that the above equation is not equivalent to 
the Euler-Lagrange equations, i.e. the constraint cannot be recovered from 
the Noether current. Notice also that the Noether currents ( |5.8| ) and ( |5.9| ) 
are different due to the difference in the corresponding Lagrangians. 



5.2 Time Translation Invar iance 



Lagrangian densities (|2.3|) and (|4.3|) are equivariant with respect to the group 
R action on Y, given by ry : (x, t, y) i— > (x, t + r, y) for any tGI. This is 
because the Lagrangians are explicitly time independent. One can readily 
compute the jet prolongation of the corresponding infinitesimal generator 
vector field Qy = (0, £, 0), where r = exp£. Then, the pull-back by j l 4> of 
the covariant momentum map corresponding to this symmetry, which we 
denote by to distinguish it from expressions in the previous section, is 
given by the following n-form on X: 



- C (e(jV)^o + Pa j (.j 1 ^ a d n x j ) (iV), 

where, in the last equality, we have used the definition of the energy density 
e given by (|2.5|). 



Noether's Theorem 5.1 implies that the exterior derivative of this expres- 
sion will be zero along solutions of the Euler-Lagrange equations. Comput- 
ing this divergence for an arbitrary £ recovers the energy continuity equation. 
For a barotropic fluid, it is given by 

' 'dcP 
dx 



e = -y/det[G\ DIV ^PJ 
while for standard elasticity the equation has the form: 




e = v / det[G] BW(V a j 4> a ). 

The expressions for an incompressible fluid and elastic medium are similar. 

Alternatively, one can consider the inverse statement and require that 
d (Jc(C)) = 0. This forces the energy continuity equation to be satisfied for 
some arbitrary section <fi. 



35 



6 Concluding Remarks and Future Directions 



In the last section of our paper we would like to comment on the work 
in progress and point out general future directions of the multisymplectic 
program. Some of the aspects discussed here are analyzed in detail in our 
companion paper Marsden, Pekarsky, Shkoller, and West [2000]. 



Other Models of Continuum Mechanics. The formalism set up in 
this paper naturally includes other models of three-dimensional linear and 
non-linear elasticity and fluid dynamics, as well as rod and shell models. For 
elasticity, the choice of the stored energy W determines a particular model 
with the corresponding Euler-Lagrange equation given by ( 2.13| ); this is a 
PDE to be solved for the deformation field cf>. Introducing the first Piola- 
Kirchhoff stress tensor V, the same equation can be written in a compact 
fully covariant form ( 2.21 ). An explicit form of the Euler-Lagrange equations 
and conservation laws for rod and shell models are not included in this 
paper but can be easily derived by following the steps outlined above. The 
constrained director models which are common in such models are handled 
well by the formulation of constraints that we use in §0. 



Constrained Multisymplectic Theories. The issue of holonomic vs. 
non-holonomic constraints in classical mechanics has a long history in the 
literature. Though there are still many open questions, the subject of linear 
and affine non-holonomic constraints is relatively well-understood (see, e.g. 
Bloch, Krishnaprasad, Marsden, and Murray [1996]). We already mentioned 
in §|| that this topic is wide open for multisymplectic field theories, partly 
due to the fact that there is simply no well-defined notion of a non-holonomic 
constraint for such theories - it appears that one needs to distinguish be- 
tween time and space partial derivatives. 

As all of the examples under present consideration are non-relativistic 
and do not have constraints involving time derivatives, we used the restric- 
tion of Hamilton's principle to the space of allowed configurations to derive 
the equations of motion. Note that this reduces to vakonomic mechanics 
in the case of an ODE system with non-holonomic constraints, and is thus 
incorrect. Of course a multisymplectic approach to non-holonomic field the- 
ories (such as one elastic body rolling, while deforming, on another, such 
as a real automobile tire on pavement), would be of considerable interest to 
develop. 
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Multisymplectic Form Formula and Conservation Laws. A very 
important aspect of any multisymplectic field theory is the existence of the 
multisymplectic form formula (2.11) which is the covariant analogue of the 
fact that the flow of conservative systems consists of symplectic maps. We 
deliberately avoid here any detailed analysis of the implications of this for- 
mula to the multisymplectic continuum mechanics and refer the reader to 
Marsden, Pekarsky, Shkoller, and West [2000], where it is treated in the 
context of Euclidean spaces and discretization. Preliminary results indi- 
cate, however, that applications of the multisymplectic form formula not 
only can be linked to some known principles in elasticity (such as the Betti 
reciprocity principle), but also can produce some new interesting relations 
which depend on the space-time direction of the first variations V, W in 
( |2.11 ). An accurate and consistent discretization of the model then results 



in so called multisymplectic integrators which preserve the discrete ana- 
logues of the multisymplectic form and the conservation laws. 



Discretization. This is another very interesting and important part of 
our project which is addressed in detail in our companion paper Marsden, 
Pekarsky, Shkoller, and West [2000], where the approach of finite elements 
for models in Euclidean spaces is adopted. It is shown that the finite ele- 
ment method for static elasticity is a multisymplectic integrator. Moreover, 
based on the result in Marsden and West [2000], it is shown that the finite 
elements time-stepping with the Newmark algorithm is a multisymplectic 
discretization. 

As we mentioned in the previous paragraph, a consistent discretization 
based on the variational principle would preserve the discrete multisymplec- 
tic form formula together with the discrete multi-momentum maps corre- 
sponding to the symmetries of a particular system. Then, the integral form 
of the discrete Noether's Theorem implies that a sum of the values of the 
discrete momentum map over some set of nodes is zero. One implication of 
this statement for incompressible fluid dynamics is a discrete version of the 
vorticity preservation. Such discrete conservations are among the hot topics 
of the ongoing research. 



Symmetry Reduction. In the previous section we discussed at length 
the particle relabeling symmetry of ideal homogeneous hydrodynamics and 
its multisymplectic realization. Reduction by this symmetry takes us from 
the Lagrangian description in terms of material positions and velocities to 
the Eulerian description in terms of spatial velocities. In the compressible 
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case one only reduces by the subgroup of the particle relabeling group that 
leaves the stored energy function invariant; for example, if the stored energy 
function depends on the deformation only through the density and entropy, 
then this means that one introduces them as dynamic fields in the reduction 
process, as in Euler-Poincare theory (see Holm, Marsden and Ratiu [1998].) 

In the unconstrained (i.e., defined on the extended jet bundle J 1 ^) mul- 
tisymplectic description of ideal incompressible fluids, the multisymplectic 
reduced space is realized as a fiber bundle T over X whose fiber coordinates 
include the Eulerian velocity u and some extra field corresponding to com- 
pressibility. Then, the reduced Lagrangian density determines, by means 
of a constrained variational principle, the Euler-Lagrange equations which 
give the evolution of the spatial velocity field u{x) £ T x together with a 
condition of u being divergence- free. A general Euler-Poincare type theo- 
rem relates this equation with equation ( [4.11 ) by relating the corresponding 
variational principles. 

Such a description is a particular example of a general procedure of 
multisymplectic reduction. The case of a finite-dimensional vertical group 
action was first considered in Castrillon-Lopez, Ratiu and Shkoller [2000]. 
More general cases of an infinite-dimensional group action such as that for 
incompressible ideal hydrodynamics, electro-magnetic fields and symmetries 
in complex fluids is planned for a future publication. The reader is also 
referred to a related work by Fernandez, Garcia, and Rodrigo [1999]. 



Vortex Methods. One of our ultimate objectives is to further develop, 
using the multisymplectic approach, some methods and techniques which 
were derived in the infinite-dimensional framework and which proved to be 
very useful. One of them is the vortex blob method developed by Chorin 
[1973], which recently has been linked to the so-called averaged Euler equa- 
tions of ideal fluid (see Oliver and Shkoller [1999]). 



Higher Order Theories. Constraints involving higher than first-order 
derivatives are beyond the current exposition and should be treated in the 
context of higher-order multisymplectic field theories defined on J k Y, k > 1. 

The averaged Euler equations (see Holm, Marsden and Ratiu [1998] and 
Marsden, Ratiu and Shkoller [2000] and references therein) provide an in- 
teresting example of a higher order fluid theory with constraints (depending 
only on first derivatives of the field) to which the multisymplectic meth- 
ods can presumably be applied by using the techniques of Kouranbaeva and 
Shkoller [2000]. It would be interesting to carry this out in detail. In the 
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long run, this promises to be an important computational model, so that 
its formulation as a multisymplectic field theory and the multisymplectic 
discretization of this theory is of considerable interest. 

Covariant Hamiltonian description. Finally, another very interesting 
aspect of the project is developing the multi-Hamiltonian description of con- 
tinuum mechanics along the lines outlined in Marsden and Shkoller [1999]. 

Acknowledgments 

The authors would like to thank Sanjay Lall, Tudor Ratiu and Michael Ortiz 
for their comments and interest in this work. 

References 

Arnold, V.I. [1966] Sur la geometrie differentielle des groupes de Lie de 
dimenson infinie et ses applications a l'hydrodynamique des fluids par- 
faits. Ann. Inst. Fourier, Grenoble 16, 319-361. 

Arnold, V.I. and B. Khesin [1998] Topological Methods in Hydrodynamics. 
Appl. Math. Sciences 125, Springer- Verlag. 

Bloch, A. M. and P.E. Crouch, [1999] Optimal control, optimization, and 
analytical mechanics. Mathematical control theory, 268-321, Springer, 
New York, 1999. 

Bloch A.M. and P.S. Krishnaprasad and J.E. Marsden and R.M. Murray 
[1996] Nonholonomic mechanical systems with symmetry, Arch. Rat. 
Mech. An., 136, 21-99 

Castrillon-Lopez, M., T.S. Ratiu and S. Shkoller [2000] Reduction in Prin- 
cipal Fiber Bundles: Covariant Euler-Poincare Equations, Proc. Am. 
Math. Soc, to appear. 

Chorin, A. [1973] Numerical study of slightly viscous flow, J. Fluid Mech. 
57, 785-796. 

Ebin, D.G. and J.E. Marsden [1970] Groups of diffeomorphisms and the 
motion of an incompressible fluid, Ann. Math. 92, 102-163. 

Fernandez A., P.L. Garcia, and C. Rodrigo [1999] Stress-energy-momentum 
tensors in higher order variational calculus, J. Geom. Phys., to appear. 



39 



Gotay, M., J. Isenberg, and J.E. Marsden [1997] Momentum Maps and the 
Hamiltonian Structure of Classical Relativistic Field Theories, I and 
II, available from http://www.cds.caltech.edu/~marsden/. 



Gotay, M.J. and J.E. Marsden [1992] Stress-energy-momentum tensors and 
the Belifante-Resenfeld formula. Cont. Math. AMS. 132, 367-392. 

Holm, D.D., J. E. Marsden and T. S. Ratiu [1998] The Euler-Poincare 
equations and semidirect products with applications to continuum the- 
ories, Adv. in Math., 137, 1-81. 

Jalnapurkar S.M. [1994] Modeling of Constrained Systems, 
http : //www, cds . caltech. edu/""{}smj/| . 

Kijowski, J. and Magil, G [1992] Relativistic elastomechanics as a La- 
grangian field theory. J. Geo. Phys. 9, 207-223. 

Kijowski, J. and W. Tulczyjew [1979] A Symplectic Framework for Field 
Theories. Springer Lect. Notes in Physics 107. 

Kouranbaeva S. and S. Shkoller [2000] A variational approach to second- 
order multisymplectic field theory, J. Geom. Phys., to appear. 

Lewis A.D. [1996] The geometry of the Gibbs-Appell equations and Gauss's 
Principle of Least Constraint, Rep. on Math. Phys., 38, 11-28. 

Lewis A.D. and R.M. Murray [1995] Variational principles for constrained 
systems: theory and experiment, The Intern. J. of Nonlinear Mech., 
30, 793-815. 

Lu J. and P. Papadopoulos [1999] A covariant constitutive approach to 
finite plasticity, to appear in Proc. 8th International Symposium on 
Plasticity. 

D. G. Luenberger [1969] Optimization by Vector Space Methods. John 
Wiley, New York. 

Marsden, J.E. [1988], The Hamiltonian formulation of classical field theory, 
Cont. Math. AMS 71, 221-235. 

Marsden, J.E. [1992], Lectures on Mechanics London Mathematical Society 
Lecture note series, 174, Cambridge University Press. 

Marsden, J.E. and T.J. R. Hughes [1983] Mathematical Foundations of Elas- 
ticity. Prentice Hall, reprinted by Dover Publications, N.Y., 1994. 



40 



Marsden, J. E., G. W. Patrick, and S. Shkoller [1998] Multisymplectic Ge- 
ometry, Variational Integrators, and Nonlinear PDEs, Comm. Math. 
Phys. 199, 351-395. 

Marsden J. E., S. Pekarsky, S. Shkoller and M. West [2000] Multisymplectic 
continuum mechanics in Euclidean spaces, preprint. 

Marsden, J.E., T.S. Ratiu, and S. Shkoller [2000] The geometry and anal- 
ysis of the averaged Euler equations and a new diffeomorphism group, 
Geom. and Fund. An. (to appear). 

Marsden, J.E. and T.S. Ratiu [1999] Introduction to Mechanics and Sym- 
metry. Texts in Applied Mathematics, 17, Springer- Verlag, 1994. Sec- 
ond Edition, 1999. 

Marsden J. E. and S. Shkoller [1999] Multisymplectic geometry, covariant 
Hamiltonians, and water waves Math. Proc. Camb. Phil. Soc, 125, 
553-575. 

Oliver M. and S. Shkoller [1999] The vortex blob method as a second-grade 
non-Newtonian fluid, E-print, 
http : //xyz . lanl . gov/abs/math . AP/99100887 . 

Shkoller, S. [1998] Geometry and curvature of diffeomorphism groups with 
H 1 metric and mean hydrodynamics, J. Func. Anal. 160, 337-365. 



41 



